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Abstract 

SU(2) Yang-Mills theory at finite extension or, equivalently, at finite temper- 
ature is probed by a homogeneous chromomagnetic field. We use a recent 
modified axial gauge formulation which has the novel feature of respecting 
the center symmetry in perturbation theory. The characteristic properties of 
the ^-symmetric phase, an extension-dependent mass term and antiperiodic 
boundary conditions, provide stabilization against magnetic field formation for 
sufficiently small extension or high temperature. In an extension of this inves- 
tigation to the deconfined phase with broken center symmetry, the combined 
constraints of thermodynamic and magnetic stability are shown to yield many 
of the high temperature properties of lattice SU(2) gauge theory. 



1 Permanent Address: Institute of Theoretical and Experimental Physics, Moscow 117259, Rus- 
sian Federation 



1 Introduction 



Although no direct phenomenological data is available for hadronic physics at high 
temperature, Monte Carlo lattice simulations of pure Yang-Mills theory and QCD 
begin to unfold a consistent picture. In pure Yang-Mills theory, a phase transition 
occurs where the Polyakov loop as order parameter signals the change from a confined 
to a deconfined phase. This corresponds to a change from the phase with unbroken 
center symmetry to the spontaneously broken phase. Moreover, there is now ample 
evidence fl], |2j of nontrivial gluonic structure above T c which persists over a significant 
temperature range. The approach to the high temperature limit of a free gluon gas 
is rather slow, and even at very high temperatures standard perturbation theory has 
only met with limited success. 

In this paper, we seek to gain some additional insight into this behaviour in pure 
SU(2) gauge theory by exploiting those nonperturbative features which are analyti- 
cally accessible in the modified axial gauge. Recently, an effective action was derived 
by some of us || for the situation in which one of the three space dimensions is com- 
pact (at "finite extension"). As a result of integrating out the gauge field variables 
which describe the phase of the Polyakov loop, this effective action has two novel fea- 
tures: a mass term and antiperiodic boundary conditions in the finite spatial direction 
for charged gluons. The mass and changed boundary conditions are themselves non- 
perturbative and a consequence of the nontrivial Jacobian in this gauge. In the same 
work, the correspondence between this finite extension formulation and the analogue 
system in thermodynamic equilibrium at finite temperature was detailed. Thus in the 
present work we shall use the terms "temperature" and "extension" interchangeably, 
with the following replacements understood, T «-> 1/L for temperature and inverse 
extension, and, respectively, e <-> — P for energy density and pressure. 

The most striking feature of the modified axial gauge is the fact that the cen- 
ter symmetry is realized in perturbation theory, implying an infinite free energy of 
a single static quark. Concomitantly, the correlator of Polyakov loops produces a 
linear potential. However, perturbation theory is clearly not adequate to describe 
the dynamics of confinement, so that for instance the value of the string tension is 
unrealistic. 

To apply this formalism to the high temperature regime is nontrivial since a 
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breakdown of the center symmetry has to occur. We cannot cross the critical point by 
using only perturbative methods. If one can at all think about the high temperature 
phase in terms of boundary conditions and masses for otherwise weakly interacting 
gluons, it is clear that the boundary conditions of charged gluons cannot remain 
antiperiodic above T c , but should, at least at very large T, become periodic. 

A useful tool for gaining some insight beyond perturbation theory is that of back- 
ground fields. One probes a given system with a specific external field and looks for 
regimes of stability respectively instability in the energy density, at a certain order 
in the fluctuations about this external field. For SU(2) Yang-Mills theory a simple, if 
crude, choice of background field is the Savvidy ansatz consisting of a constant colour 
neutral chromomagnetic field [|J. This ansatz itself has a long history at both zero 
and nonzero temperatures. Indeed, various authors have tried to exploit the energy 
density or effective potential for a Savvidy type external field in order to observe the 
deconfinement phase transition, by studying the change in the position of the mini- 
mum of the energy density. Results, however, have been under dispute to date due 
to the varying degrees of ad hoc measures used in these works. In the original (zero 
temperature) study of Savvidy, the minimum of the energy density is driven away 
from the zero field limit by a quite specific mode; the same mode however makes the 
potential imaginary at the local minimum that is finally established (to one loop), so 
that this minimum is in turn unstable. At finite temperature for the standard Yang- 
Mills action, the unstable mode persists: a number of authors ||, || have simply cut 
this mode out and obtained a nontrivial minimum at low temperature which becomes 
a minimum at zero field for high temperature. This was then interpreted as a signal 
of the deconfinement phase transition. However, as demonstrated in J7|, , even with 
quarks once the unstable mode is put back in, at high temperature the instability 
persists: there is no temperature at which a stable minimum to the one loop energy 
density at zero field develops in the background field gauge. In another approach |J, 
and one quite similar to ours, this effective potential has been studied by treating the 
boundary conditions of charged gluons as a variational parameter. For any nonperi- 
odic choice of boundary conditions, the corresponding shift in the lowest Matsubara 
frequency is found to make the theory stable at sufficiently high temperatures. How- 
ever, with the change in boundary conditions in the field dependent part only, these 
results are hard to interpret physically. 
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Our strategy will be as follows. After reviewing the key results of 0, we shall study 
the Z2 symmetric phase, that is the phase with the center symmetry realized, of the 
SU(2) Yang-Mills theory via its response to an external magnetic field of the Savvidy 
type. To this end, we shall compute the energy density in zeta-function regularization 
JTIJ for SU(2) gluons in the presence of the Savvidy background field. For later 
generalization we perform this calculation for arbitrary mass and boundary conditions 
of the charged gluons. With the "canonical" masses and boundary conditions as 
derived in the modified axial gauge theory in |§, the only minima in the energy 
density that appear are at zero external field, but for sufficiently high temperature 
these minima are stable. The Z2 symmetric, perturbative phase is however not and 
cannot be the correct high temperature phase of the SU(2) Yang-Mills theory. If 
antiperiodic boundary conditions persisted beyond the "critical" point, dimensional 
reduction would leave a 2 + 1 dimensional Maxwell theory rather than the non- 
Abelian theory that is seen in lattice computations. Furthermore, for sufficiently 
high temperatures where we do not expect significant non-perturbative dynamics, the 
thermodynamic instability of the Z 2 symmetric phase becomes an important issue. We 
therefore will study in the final part the consequences of magnetic and thermodynamic 
stability relevant for the high temperature phase. Allowing for variations in mass 
and boundary conditions, it will be possible to simultaneously satisfy both stability 
requirements. As displayed in a discussion of the high temperature properties of the 
system, the resulting phenomenological description shares, qualitatively, important 
characteristics with finite temperature lattice gauge theory. 



2 Review of the Modified Axial Gauge 

At finite extension or temperature, one space-time direction is singled out; hence, 
an axial type of gauge is particularly convenient. To be specific, we choose the 3- 
direction as the compact direction (x^ G [0, L\) and work in Minkowski space. All our 
results can equally well be reinterpreted in terms of finite temperature equilibrium 
thermodynamics if desired, by going to Euclidean space and exploiting covariance. 
The other space-time directions will be denoted by a, (3 = 0, 1,2 and referred to as 
"transverse" . As is well known, periodic boundary conditions imposed on gauge fields 
in the 3-direction do not permit complete elimination of A3; a gauge invariant zero 
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mode has to be retained. More precisely, this 3-dimensional residual field corresponds 
to the eigenphases of the field W(x±) = P exp(ig Jq dx 3 A 3 ), the trace of which can 
be identified with the Polyakov loop, w(xj_) = |Tr W(x±). We thus demand only the 
weaker, "modified" axial gauge condition d 3 A 3 = 0, which can be supplemented by 
a diagonalization in color space, A 3 —>■ a|r 3 /2. For simplicity, we denote a 3 as 03 in 
the following. The Polyakov loop is then given by 

w (x±) = cos (gLa 3 (x±) /2) , (1) 

showing that in the modified axial gauge this important variable appears as an ele- 
mentary rather than composite field. 

The presence of the three dimensional field a 3 (x±) indicates the necessity for a 
further gauge fixing. Indeed, gauge ambiguities are still present in the sector of 
transverse (a = 0,1,2), neutral, x 3 independent gauge fields. Due to their Abelian 
nature, these ambiguities can easily be resolved by adding a standard gauge fixing 
term to the Lagrangian. The 3-dimensional Feynman gauge is particularly convenient 
for perturbative calculations, 

= f L dxsPA*{x)\ . (2) 
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Having fully specified the gauge, we arrive at the generating functional 

Z [JJ 3 ] = J D [A a ] d [a a ] exp |^ YM [A a , a 3 ] + tS gi [A 3 a ] + 1 J d*x (A a a J a ' a + a 3 j 

(3) 

(see also Ref. Here, SyM(A a , a 3 ) is the standard Yang-Mills action with A 3 

replaced by a 3 ; the integration measure for the 03 functional integral is given by 

d [a 3 ] = [] sin 2 (gLa 3 (y ± )/2) 6 (a 3 ( Vl _)) (2n/gL - a 3 (yj_)) da 3 (y ± ) (4) 

y± 

and accounts for the fact that gLa 3 /2 appears in the parametrization of the group 
manifold S3 as the first polar angle. As a consequence the corresponding part of the 
Haar measure enters with a finite range of integration. In the canonical formulation, 
this Faddeev-Popov determinant arises as a Jacobian modifying the kinetic energy 
of the Polyakov loop variables a 3 (x±) fl2"fl . In QED the same procedure yields the 
standard flat measure for a 3 (x±). 
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The finite range of integration for 03 in the above functional integral precludes 
a conventional perturbative treatment, since the quadratic part of the action does 
not yield a Gaussian integral. Here, we follow the approach proposed in Ref. ||, 
which amounts to integrating out 03 completely prior to any perturbative expansion. 
First, we regulate the functional integral by a transverse lattice (lattice constant €). 
Subsequently, the variable 03 is shifted to a' 3 = as — ir/gL. As a consequence, the 
sin 2 factor in the measure (^) is replaced by cos 2 , and the integration limits become 
[—7r/gL,7r/gL]. In order to preserve the standard form of the minimal coupling of 
charged gluon fields to a' 3 , it is advantageous to redefine these fields as follows, 

A l =f iA 2 -> exp(±i7rx 3 /L)(A! =f iA 2 ) ■ (5) 

This does not change the physics, but yields new charged gluon fields which, in 
contradistinction to the neutral ones, obey antiperiodic boundary conditions. 

The fact that one can integrate out completely rests on the following observa- 
tion: as a consequence of the finite range of integrations, the field 03 (and hence the 
Polyakov loop) becomes "ultralocal" ; it does not propagate in the continuum limit 
0. In lattice terminology, one finds that hopping terms are suppressed by factors 
lj g 2 L. Since the coupling constant is expected to "run" with I in the standard loga- 
rithmic manner (as will indeed be confirmed below in the present framework) and L 
is a macroscopic length, this factor goes to zero in the continuum limit. 

Formally, ultralocality of the Polyakov loop field 03 means that it can be com- 
pletely integrated out. Introducing a lattice regularization in the transverse direc- 
tions in order to define the functional integral over this field, one arrives at finite 
non-Gaussian integrals which can be explicitly carried out. No ultraviolet divergent 
expressions occur in this subsector of the theory, so that the lattice regulator can be 
safely taken to zero, leaving an effective action for the remaining degrees of freedom 
(a = 0,1, 2): 

r 1 M 2 r 

S efi [A a ] = S YM [A a , A 3 = 0} + S gi [Al] + -f £ / d A xA a a (x) A a ' a (x) . (6) 

^ 0=1,2"' 

In this effective action, the Polyakov loop variable has left its signature in the "geo- 
metrical" mass term of the charged gluons 

Ml = (tt 2 /3 - 2) /L 2 (7) 
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and in the change to antiperiodic boundary conditions 

A 1 a 2 (x ± ,L) = -A 1 a 2 (x ± ,0) , (8) 

while the neutral gluons A z a remain massless and periodic. The antiperiodic bound- 
ary conditions are a consequence of the mean value of the Polyakov loop variable, 
the geometrical mass comes from its fluctuations. Notice that both of these changes 
in the action are independent of the coupling constant. In the center symmetry 
and a properly redefined order parameter for the confinement-deconfinement tran- 
sition in this effective theory are discussed in detail. Here we only recall that the 
expectation value of the Polyakov loop (w(x±)) vanishes in zeroeth order perturba- 
tion theory, indicating that this formulation is particularly well suited for describing 
the Z 2 -symmetric phase. 

Having reviewed the essential aspects of 0] it is useful to look at this effective 
theory in a wider context. The fact that the Polyakov loops are indicating confine- 
ment of static quarks means that we have a natural formalism for dealing with finite 
temperature gauge theories which are still in the confined phase. As the tempera- 
ture is increased, one knows from lattice calculations that a phase transition to a 
deconfined phase happens at a certain critical temperature T c ; correspondingly, the 
Polyakov loop order parameter acquires a non-vanishing expectation value. We can- 
not expect to be able to treat this transition perturbatively. However, we can try to 
extrapolate our findings into the deconfined phase. If there exists a similarly ade- 
quate kind of perturbation theory for the broken phase where gluons can be treated 
as weakly interacting, it should not involve antiperiodic boundary conditions, but 
rather quasi-periodic ones which approach periodic boundary conditions in the limit 
T — > oo. This is just another way of phrasing the fact that Polyakov loops acquire an 
expectation value in the deconfined phase. Indeed, a theory of essentially free gluons 
with antiperiodic boundary conditions for very high temperatures can be ruled out on 
the grounds of dimensional reduction: in the limit T — > oo for any choice of boundary 
conditions other than periodic, charged gluons, exhibiting a gap in their single par- 
ticle spectra, would decouple and the resulting dimensionally reduced theory would 
not be three-dimensional Yang-Mills theory but rather Maxwell theory, in conflict 
with lattice results. Another argument comes from the sign of the Casimir pressure, 
which for antiperiodic boundary conditions would signal thermodynamic instability 
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at high temperature, an unacceptable result. Finally, above the critical temperature 
the effective mass of charged gluons cannot increase forever linearly with temperature, 
both from dimensional reduction arguments and in order that the Stefan-Boltzmann 
law can be truly restored at very high temperatures. 

To summarize, nontrivial boundary conditions, a (magnetic) gluon mass and the 
realization of center symmetry at finite temperature are the basic non-perturbative 
features brought about by gauge fixing in our formalism. Now, as the temperature 
is increased a phase transition to the broken phase must occur. At this point, evi- 
dently the boundary condition angle and the gluon mass must change. The physics 
associated with this change is computationally difficult to describe in our approach, 
but nonetheless this change in the values of these parameters will not be immedi- 
ate: namely nontrivial values should persist above the phase transition. The simplest 
way then to characterize the symmetry broken phase is to use these parameters, 
though now allowing them to take values other than those determined above. The 
requirement of stability combined with the demand that the Stefan-Boltzmann law 
be approached at high temperatures can enable us to study the behaviour of these 
parameters in the very high temperature regime. We thus have variables x an d M as- 
sociated with charged gluons in the effective theory of 0. The first quantity describes 
the angle in the boundary condition 

A 1 a 2 (x L ,L) = e^A 1 a 2 (x u 0) (9) 

with < x < 27r, the second the mass of these gluons. In terms of these parameters 
we calculate the energy density for an external magnetic field of the Sawidy type. We 
shall then set these parameters to their values x = 71 1 M = M3 for charged gluons, 
X = M = for neutral gluons, which characterize the ^-symmetric phase, and 
investigate the stability problem. In the next step, we shall treat these parameters in 
a less restrictive fashion in order to explore the high temperature phase. 

3 Energy Density in a Background Magnetic Field 

In the above formalism where the redundant variables have been eliminated, we pro- 
ceed to study the stability properties of the effective theory of Eq. @ by probing it 
with a background magnetic field. We choose a homogeneous external chromomag- 
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netic field in color 3-direction, with the potential of the form 

^Ib g = -5 a3 g^x 2 H . (10) 

Although this looks identical to Savvidy's choice Q at first glance, we should like 
to point out an important difference. Traditionally, the background field has been 
introduced before gauge fixing. In this context, the choice of the background field to 
be diagonal has no physical meaning. In our case however, the 3-direction has already 
been specified by diagonalizing the Polyakov loop variables. What we investigate then 
is the stability of the vacuum against homogeneous magnetic field fluctuations which 
point in the same direction in color space as the Polyakov loop. 

The energy density in the presence of the external field is the sum over the single 
particle energies. Since we cannot use the popular background gauge, we must solve 
the Landau level problem in the axial gauge; some details of the calculation can be 
found in the appendix. Clearly, the neutral gluons contribute only to the Casimir (H 
independent) part of the energy density, which will be included later. We thus here 
restrict ourselves to the charged gluon contributions. The expressions at this point 
are ultraviolet divergent and in need of regularization which will be done here by 
zeta-function methods. In this scheme, we need to evaluate the following expression 
for the vacuum energy density, 

H 2 . </B 



I ZTTL s =±ln=0k=-<x 



2 gH(s + n + 1/2) + (2?rA: r | x)2 + M 2 1 



The first term is just the classical energy density of the external magnetic field. The 
coefficient of the second term comes from the density of states of Landau levels, 
the parameter e regularising the ultraviolet behaviour of the sum. The arbitrary 
scale parameter /i 2 has been introduced to keep the dimensions correct. The Casimir 
contribution from the charged gluons will emerge in the H — > limit of this expression. 
Eq. (|ll|) can be evaluated using standard techniques |10| the details of which we omit 
here. We shall give later the full result of the summations after renormalization of the 
ultraviolet singularities, for the purpose of which it is here sufficient to just consider 
the weak field limit of the result gH <C M 2 , 

IP llg 2 H 2 
£ ~ ~2~ + 48vr 2 ["' \ /i 2 J e ' ' ' r V 2, 



In ( % ) - i + l + ^(-i)-4£ cos(n X )K (nML) 



n=0 
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M 4 I" ( M 2 \ 1 , / 1\ 11 2M 2 ~ cos(mx) T , , „, r , , , 

We recognize in this the ultraviolet divergent part (pole in e). This itself has pieces 
dependent on H and on M, L. Also, the H independent part of this expression is the 
unrenormalized Casimir energy for two free scalar boson fields. We consider this in 
some detail first now. 



3.1 Casimir energy and pressure 

Setting H = and using that ip(—l/2) = 2 — 7 — 2 In 2, we obtain straightforwardly 

+ + H^/M)) (ML/2)" 



e(H = 0) = 2 



1 3 I 

00 

n 2 



VT 2 L 4 

+ (ML) 2 ]T cos{n X )K 2 {nML) 



n=l 



(13) 



To eliminate the divergence for e — > 0, we have to properly renormalize the energy 
density. If M ~ 1/L and x is -^-independent the energy density can be calculated 
with respect to an arbitrary reference point L, as follows 

The logarithmic term in Eq. ( |i~4"| ) prevents approach to the Stefan-Boltzmann law at 
small L and also modifies the standard relation between pressure and energy density 
(remember that e and — P are interchanged as compared to the thermal situation) 

p(l) = -JL (£L)=HL) + ^. ( 15 ) 

If on the other hand we would simply ignore the mass, but work with antiperiodic 
boundary conditions x = ^1 the result would be finite but exhibit a repulsive Casimir 
effect, signalling thermodynamic instability of the corresponding finite temperature 
system, 

P(L) = ML) = ~ • (16) 

This illustrates the above mentioned difficulties, if we would apply the "canonical" 
values of x an d M as inferred from the axial gauge in the Z2 symmetric phase also 
at small L (high T). 
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3.2 Field dependent part 

With the external field non-zero, the ultraviolet divergence contains both field de- 
pendent and M, L dependent parts. The renormalization must proceed in two steps: 
first a coupling constant renormalization after identification of the renormalization 
group invariant B = gH is performed. This leads to the correct SU(2) Yang-Mills 
beta function to one loop order, a confirmation that the non-perturbative aspects of 
the present formulation do not spoil asymptotic freedom. The arbitrary scale \j? is 
fixed by choosing 
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2gl(fi) 48vr 



(17) 



where Ams is the scale parameter for SU(2) in the minimal subtraction scheme. There 
remain now poles independent of B, but dependent on M, L. These are eliminated 
by performing a subtraction at B = 0. Finally, for arbitrary field B we obtain the 
rather lengthy renormalized expression: 



e K (L,B) = e{L,B)-£{L,0) 



+ 



+ 



+ 



11B 2 
48tt 2 



In 



2B 
A 2 



B 1 



8tt 2 



In ( 





2B 






167T 2 \ 


M 2 


M 2 + B 


)* 


2\B\ 





2tt 2 
2B x 



M 2 + B" 
2B 



+ V>(-l/2)l+(fl- 



-B) 



J2[B(2n + 1) + M 2 } 1 ' 2 J <^pd Kl ( mL [ B (2n + 1) + M 2 ] 1/2 



7l 2 L 

B(M' 2 + B) 



n=0 

2 i dU/2 oo 



m=l 



m 



n 2 L 



X; ^^-K^nLiM 2 + B) 1 ' 2 ) + (B^ -B) 



n=l 



n 



M 4 . , . . 2M 2 ~ cos(mv) . 

(V(-l/2) - 1/2) + — X -^lK 2 (mML) 



16n 



7T 



2 L 2 



m=l 



77T 



18) 



with the Riemann zeta-function £(s,i/), P~5[] , and £'(s,v) = d£(s,v)/ds. This result 
is more involved than that of previous authors mainly due to the non- vanishing gluon 
mass. 

We consider two limits of this renormalized expression, the weak field (B <C M 2 ) 
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or the strong field case (B ^> M 2 ). Taking the strong field limit first, we obtain 

11B\ ( B 

L MS/ 

with subleading terms having an imaginary part now. This is precisely the case 
considered by Savvidy and many others. In the weak field case, we obtain 



in | (19) 



1 1 R 2 



48vT 2 



( M 2 \ 00 
In — - + 1 + ^(-1/2) - 4 J2 cos(n X )K (nML) 



A 2 , 

IV MS/ n=l 



(20) 



Note that due to the existence of the additional mass scale, we no longer have the 
magnetic field appearing in the logarithm. This has important consequences for the 
discussion of magnetic stability. 

We now input the "theoretical" values characteristic for the Z 2 symmetric phase 
for the parameters M and x-, i- e - M = M 3 = (7r 2 /3 — 2) 1,/2 /L, \ = 77 (note that only 
the charged gluon contributions have survived after the zero-field subtraction). We 
have numerically computed £r for various values of the extension L against B. As 
a result, we find no sign of minima at nonzero i?-field values for which there is no 
imaginary part. Only at the point B = can a stable minimum occur. Comparing 
the numerical to the weak field result demonstrates the effectiveness of the lowest 
order B contribution in the neighborhood of the minimum. This is shown in Fig. 1. 
For larger values of L, the energy density would start with a negative curvature and 
the cusp visible in Fig. 1 be more pronounced. At the cusp, the energy develops an 
imaginary part; this will be discussed further in the following subsection. 

Having established that B = represents the only relevant minimum, it is simple 
to determine the critical value of L below which the sign of B 2 in Eq. (|20|) is positive 
for the canonical values for \ and M. It is found to be 

L c = —^exp[(l + ^(-l/2) -Aj2(-l) n K (nLM 3 ))/2] « ^ . (21) 
Ams n =i a ms 

For magnetic stability, one must have L < L c . For these values of L the system can be 
said to be preferring the "empty" vacuum (B = 0), even though the calculation does 
not yet contain enough dynamics to make statements about the true small L (high 
temperature) behaviour. Given that T c w Ams> this means that the center-symmetric 
phase is stable above roughly T c /3. Thus, anticipating the final section, it is plausible 
to assume that in the broken phase this stability will persist. 
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3.3 Imaginary part of energy density 

Above, we have focussed on the sign change of the B 2 contribution to the renormal- 
ized energy density, i.e., the weak field aspect. Our calculation shares with previous 
investigations the disease that above a certain critical magnetic field, the energy den- 
sity develops an imaginary part. In fact, this happens at all values of L. In Fig. 1, it 
is reflected in the cusp in the real part of e, whose position corresponds to the thresh- 
old magnetic field; the curves at smaller L develop such cusps at larger I?-values not 
visible in the figure. The value of the critical magnetic field can easily be understood: 
this is nothing but the point at which the mass gap (due to the mass and nonzero 
lowest Matsubara mode) is completely compensated for by the magnetic field (the 
lowest Landau level of the s = — 1 gluon reaches zero, see Eq. (|TTD). This kind of in- 
stability at large magnetic fields is more problematic than the one at weak fields, since 
we cannot trust the one loop calculation anymore. Clearly, at the critical magnetic 
field, a restructuring of the vacuum must take place since the gluons in the lowest 
Landau level can condense. Here, we have made no attempt to treat these particular 
modes non-perturbatively and therefore cannot handle this kind of instability. One 
interesting aspect however is worth mentioning, which is specific for our approach: 
as elaborated in ref. M, the Polyakov loop correlator is dominated at large distance 
by the closest singularity in momentum space, i.e. the threshold for producing two 
charged gluons. On the other hand, it is related to the potential between static fun- 
damental sources. On the basis of the Landau orbits, we can immediately predict 
that the "string tension" will go like 



with c 2 = 47r 2 /3 — 2 coming from the combination of the mass term and the lowest 
Matsubara mode. At the cusp, the string tension vanishes; this opens the possibility 
for exploring the center symmetry breaking as a function of a parameter other than 
the temperature, namely the strength of the background field. We refrain from doing 
this here since it would require a better dynamical treatment of those modes which 
might condense. 




(22) 
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4 High Temperature Regime 



We appear to have arrived at a partly satisfactory description of the response of the 
Z 2 symmetric Yang-Mills phase to a homogeneous magnetic field at some interme- 
diate temperature. As emphasized above, extension towards high temperatures is 
problematic. Although increasing the temperature improves the stability with re- 
spect to magnetic field fluctuations, the thermodynamic instability now becomes an 
important issue. Criteria of thermodynamic stability generally imply second order 
derivatives of thermodynamic quantities. We expect at sufficiently high temperature 
a gas of gluons in a regime where the Stefan-Boltzmann law is approximately sat- 
isfied, namely the T 4 dependence of energy density or pressure is dominant. Thus 
our simplified criterion for stability will be positive pressure for thermodynamics, or 
equivalently negative energy density for finite extension. 

Now, the requirements of stability against magnetic field formation and ther- 
modynamic stability are to some extent in conflict. Roughly speaking, magnetic 
(thermodynamic) stability demands large (small) values of the "infrared parameters" 
M, x- It is thus of interest to discuss the consequences of the combined constraints 
of magnetic and thermodynamic stability. This study will be performed by consider- 
ing variations in the quantities M and x i n order to achieve simultaneously stability 
with respect to both thermal and magnetic field fluctuations. To allow for variations 
in these parameters is clearly appropriate for a description of the high temperature 
phase. At temperatures or extensions beyond the confinement-deconfinement transi- 
tion with the center symmetry spontaneously broken, the quantities M and x must 
deviate from their canonical values characteristic for the Z 2 symmetric phase. The 
aim of our studies is to display the consequences of the stability rather than to at- 
tempt a fit to lattice results for the high temperature phase. We will not allow for 
simultaneous parameter changes, but identify one of the two parameters with its infi- 
nite temperature limit (M = or x — 0). In the first place, we shall assume massless 
charged gluons with L-dependent angle x (quasi-periodic gluons); then, we shall keep 
the boundary conditions of the charged gluons periodic but allow for a more general 
L-dependence of their mass (massive gluons). 
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4.1 Stability with quasi-periodic gluons 

Here we assume ML to be negligible. We first consider the requirement of thermo- 
dynamic stability. The Casimir energy density for the massless charged gluons with 
quasi-periodic boundary conditions can be deduced from Eq. (|13|) . 

4 ^ cos(nx) 4tt 2 / x x 



vr 2 L 4 ^ n 4 3L 4 * V2tt 

with the Bernoulli polynomial B&(x) = — 1/30 + x 2 (l — x) 2 . Eq. (p3|) displays the 
Stefan-Boltzmann behaviour for periodic boundary conditions as well as a change of 
sign for antiperiodic ones. The Casimir pressure is given by 

Positive pressure in thermodynamics corresponds to a negative Casimir energy den- 
sity, Eq. (|23"D , and therefore requires 



X < 1-51 • (25) 

A positive susceptibility is a necessary condition for magnetic stability. The vanishing 
of the B 2 term in the energy density for the case ML = can be determined by 
carrying out the sum in Eq. (|20| ) with the help of the identity 

T,™s(n X )K (ny) = - (7 + In + ^ + ^ + (26) 



7T ~ / ll 1_ 

V v /y 2 + (2t™-x) 2 yV + (27rn + x) 2 7rn / 
The requirement that expression (|20|) vanishes then simplifies to 

LAMs = 4vrexp|^l-^ + ^(-l/2)+^(l + ^)+^(l-^))} - (27) 

This relation, together with the inequality (^), divides the {x,L) plane into mag- 
netically stable and unstable regions as seen in Fig. 2. Notice that the combined 
constraints from magnetic and thermodynamic stability are quite restrictive. Invert- 
ing Eq. (|27|) yields a lower bound for x as function of L. At small L in particular, 
this inversion can be performed in closed form with the result 

^ L ~* °) w ^71 ^ 9 w 9 7777 ^ • (28) 

3(1 - 7) + 2 ln(27r/(LA M s)) 
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It is instructive to reinterpret this formula as high T behaviour and eliminate the 
logarithm in favor of the (one loop) running coupling constant, 

<T 2 (T) = -^ln(27rT/A MS ) ■ ( 29 ) 
Then, for large T we find the magnetic stability bound 

X (T) > ^-g\T) , (30) 

now in a more conventional, seemingly perturbative, guise. In the case of quasi- 
periodic gluons, any angle \ vanishing faster then g 2 {T) (that is, faster than 1/ ln(T)) 
for T — > oo would be ruled out. 

In the absence of a complete theory for x{L), we simply identify with its 

bound obtained from the requirement of magnetic stability (Eq. ( p7|) and, at large 
T, Eq. (0)). If we then reinterpret L as 1/T and interchange e with — P (and vice 
versa) , we can get an upper limit for pressure and energy density of a gluon gas at high 
temperature, at least within the limited set of configurations which we can describe. 
This has been done in Fig. 3, where e and P are normalized to the Stefan-Boltzmann 
values e = 3P = 7r 2 T 4 /5 and the contribution from neutral gluons (which we assume 
to obey periodic boundary conditions) has been added. 

We first note that we cannot perform such a calculation at temperatures below 
T = 0.34A M s, simply because Eq. fl27p ceases to have a solution (the right hand 
side is bounded from above). At low temperatures in Fig. 3, we observe the onset 
of thermodynamic instability where the pressure changes sign (or the Casimir effect 
changes from attraction to repulsion) because x crosses the value 1.51, see Eq. (|25| ). 
The most conspicuous feature of Fig. 3 is however the slow, logarithmic approach 
to the Stefan-Boltzmann limit, reminiscent of the lattice data 0] but quantitatively 
farther away from this limit in the covered temperature region. A rather sensitive 
measure for interaction effects is the quantity e — 3P, which is of interest also in 
view of its direct relation to the trace of the energy-momentum tensor and the gluon 
condensate @. Although our predicted temperature dependence is too rapid, this 
difference has the right order of magnitude as compared to the lattice data of as 
illustrated in Fig. 4 (curve labelled "x")> we nave assumed A MS m T c lattlC0 290 MeV 
for this particular comparison. We finally note the asymptotic behaviour of the 
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thermodynamic functions shown in Fig. 3, 

£sb(T) P sb (T) 3 I A MS / ^ ' 

The difference 5 — 3P on the other hand behaves like T 4 (ln(27rT/AMs))" 3 for large 
T. 

4.2 Stability with massive gluons 

In this section, we shall set x = f° r &U L an d study the effect of an L-dependent, 
"magnetic" gluon mass. The condition for the critical mass above which the system 
is magnetically stable reads (see Eqs. ( |20| , p6|) ) 



21 n ^ = l-2 7 + «- 1 /2)-|L-2.£|^_^_--L] . (32) 



\y/(ML) 2 + (2vm) 2 ™ t 



The (numerically obtained) boundary between stable and unstable regions in the 
(M, L)-plane for this massive gluon case is shown in Fig. 5. Stability is only possible 
for values of M larger than a certain limiting value, 

Mo = A M se)q>(-i(l + V(-l/2))) ~ 0.60A MS • (33) 
The small L behaviour of M can again be derived analytically, 

^- 0|B Wi-,| t ;NW (34) 

Inserting as above the running coupling constant and interchanging L and 1/T, this 
corresponds to the lower bound 

M(T) > ^Tg\T) (T - oo) (35) 

at high temperatures. It is remarkable that this coincides with both theoretical ex- 
pectations and lattice results for the temperature dependence of the gluon magnetic 
mass; thus for instance, a recent determination in the Landau gauge in a wide tem- 



perature range |L5| has been fitted with the formula 

M m (T) = (0.46 ± 0.01) Tg 2 {T) , (36) 
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surprisingly close to (but larger than) our lower bound from magnetic stability, Eq. 
(|35|). We emphasize that we have obtained Eq. (^) using perturbative methods, 
although the magnetic gluon mass is well known to be a genuinely non-perturbative 
quantity. This is obviously a result from not calculating M directly, but inferring it 
indirectly from magnetic stability considerations. 

We can again estimate thermodynamic functions by inserting the critical mass into 
the expressions for the Casimir effect above, Eq. ([13]). The important feature here is 
that the gluon mass rises more slowly than linearly with temperature. This is highly 
welcome, since it guarantees that we will asymptotically recover Stefan-Boltzmann 
behaviour of an ideal gluon gas. Since on the other hand, magnetic stability does 
not permit a mass growing more slowly than with T/ln(T), these restrictions leave 
very little room for a different behaviour. In contrast to the massless case with In- 
dependent boundary conditions (quasi-periodic gluons), we now need one subtraction 
to get a finite result, see Eq. ([13]). This reduces somewhat the predictive power of 
the approach based on massive gluons. Nevertheless, encouraged by the reasonable 
T-dependence of the gluon mass, we proceed, adjusting the energy density at one 
reference length L. After this subtraction, we have to evaluate the Casimir energy 
density for charged gluons 

= M*(±- - 1 — 2 V 1 ( K ^ nML ) K 2 {nML) \\ 
£ U 4 16vr 2ll M K 2 hn 2 \ (ML) 2 (ML) 2 ))' 1 J 

where M = M(L),e = e(L). As a sample calculation, we set e equal to a certain 
fraction of the Stefan-Boltzmann value (0.8) at the highest temperature; this value 
has been selected since it is consistent with the reduction of 30% at 2T C reported in 
the SU(2) lattice calculation of fl6 |. The pressure is evaluated by numerical differen- 
tiation, see Eq. (|T~5D . The result is shown in Fig. 6 where a massless neutral gluon 
contribution has been added, since our formalism gives no hint to modifications in the 
neutral sector. While different in detail, this particular calculation leads qualitatively 
to a picture similar to the one with quasi-periodic gluons shown in Fig. 3. Once 
again we observe the slow (logarithmic) approach to the Stefan-Boltzmann limit as 
T —>■ oo. We have also calculated the interaction measure e — 3P in the region where 
lattice data || are available and included it into Fig. 4 (curve labelled "M" ) to test 
the sensitivity of these calculations to details. The same renormalization condition 
as in Fig. 6 has been used here. 
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4.3 Electric screening mass 

It is tempting to use the function x{T) derived from the condition of magnetic stability 
also to improve other perturbative predictions. One key quantity in hot gauge theories 
is the Debye screening mass, related to the zero momentum limit of the vacuum 
polarization tensor U 00 . Standard (finite extension) perturbation theory in the axial 
gauge gives for n 33 the lowest order expression 

-e! = -^E((2™ + x) 2 + (ML) 2 ) 1/2 , (38) 

71 -Li n 

where the gluon loop and tadpole are included and dimensional regularization has 
been used to perform the 3-dimensional integrations. In the limit M — > and evalu- 
ating the divergent sum with zeta function methods, we find 

< = %) (39) 

(B 2 (x) = 1/6 — x(l — x)). Incidentally, this result can be checked easily by taking 
the second derivative of e, Eq. (^), with respect to x- With the special case of 
periodic boundary conditions and replacing L by 1/T, we recover the well known 
result = 2g 2 T 2 /3, where g 2 should be replaced by the running coupling constant. 
On the other hand, insertion of the critical value for x(T), the x-dependent factor in 
Eq. fl39| ) changes the temperature dependence of m e i significantly. Over a wide range 
of temperatures, we observe a compensation between the decrease due to the run- 
ning coupling constant and the increase due to the restoration of periodic boundary 
conditions; see Fig. 7 (curve labelled "x"). The resulting Debye mass has a linear 
T-dependence, m e i ~ 0.84T, to very good accuracy from T = 4Ams up to T much 
larger than 16Ams shown in the plot. This temperature dependence is consistent 
with the flat behaviour observed in the lattice calculation [O] , although the prefactor 
there was found to be significantly larger than ours, m e \ ~ 2.5T. Very recently, a 
new lattice calculation |T7| reported somewhat smaller values of m e i in this tempera- 
ture region. However, for extremely high temperatures (up to 10 T c ), the logarithmic 
running was observed, although with a coefficient significantly larger than the naive 
perturbative expectation. 

In the other case where the magnetic gluon mass M serves to stabilize the system, 
we can go through similar considerations. Using zeta-function regularization, m 2 j 
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in the massive gluon case can be treated along similar lines as the Casimir energy 
density. A subtraction of m^/M 2 at a certain L is necessary to get a finite answer. 
We find (in analogy to Eq. (|37|)) 



\M 2 vr 2 M M^ in \ ML ML ) ) 1 1 

If we now adjust m e \ to a reasonable value, we obtain a behaviour similar to that 
for quasi-periodic gluons, namely a striking plateau in m e \/T up to high T (Fig. 7, 
curve labelled "M"). On the other hand, it does not seem possible to remove the 
discrepancy in the absolute value of a factor of 2 with the lattice data Jl5|] by simply 
adjusting our electric mass to a higher value at one temperature; if we were to do 
this, the curve for m e \ would approach the asymptotic limit from above, and the flat 
behaviour of m e \/T would be destroyed. Thus we cannot reconcile absolute value and 
T-dependence of the Debye screening mass better than to the shown qualitative level. 
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■■ 2 g 2 M 4g 2 ~ 1 ( K x {nML) K x (nML)\ 



el 



5 Conclusions 

In summary, we have probed SU(2) Yang-Mills theory at finite extension by means of a 
homogeneous magnetic field. Technically, our study differs in two important respects 
from previous works devoted to this topic: On the technical side, we have worked 
in a completely gauge fixed framework from the outset, the modified axial gauge. 
The gauge fixed theory exhibits certain non-perturbative ingredients like a mass and 
boundary condition angle for charged gluons in the Z 2 -symmetric (confined) phase. 
^From the physics point of view, our investigations focussed then on the interplay 
between magnetic and thermodynamic stability. The calculation then proceeded in a 
straightforward way, by evaluation of the energy density to one loop order (we prefer 
the Casimir ground state to a system in thermodynamical equilibrium, but the two 
are equivalent). A tachyonic instability at low temperature against large fields in- 
deed appeared in our calculation, But because this is in a regime where the one-loop 
calculation cannot anyway be trusted, this was not relevant for the analysis we subse- 
quently carried out. Indeed, unlike in standard background gauge computations @, U, 
in our gauge-fixed formulation of the Z 2 symmetric phase with antiperiodic bound- 
ary conditions and the mass ~ 1/L for charged gluons, we obtained stability against 
weak field magnetic fluctuations below a certain extension or, equivalently, above a 
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certain temperature. This temperature was found to be approximately one-third the 
confinement-deconfinement transition temperature T c as determined by lattice gauge 
calculations. 

The fact that our calculation achieved stability already before the actual transition 
to the Z 2 -broken phase gave us some confidence in extrapolating our formalism into 
this broken phase. To use our formulation there we made one assumption: that 
above the phase transition the fundamental nonperturbative quantities remain the 
mass and boundary condition angle of the charged gluons. In the course of the phase 
transition, their values should, naturally, change and of course we cannot predict how 
they should change. However a complete change to periodic boundary conditions 
and massless gluons at the critical point is not possible because it would lead to 
magnetic instability. On the other hand, at sufficiently high temperatures our one 
loop computation, with these minimal nonperturbative modifications, should become 
more reliable. We have determined the value of the mass, or the boundary condition 
angle, which satisfies in a minimal way the condition of magnetic stability. These 
bounds were in turn reinserted into various thermodynamic quantities. We obtained 
features which are stongly reminiscent of results of lattice calculations, but which 
have thus far resisted any kind of perturbative understanding. Specifically, we found 
that the slow approach of pressure and energy density to the Stefan-Boltzmann limit, 
seen in lattice results [H , within our treatment seems quite unavoidable. Even more 
surprising is our success in reproducing the linear T— rather than gT— dependence for 
the electric screening mass seen by \\L5\ . We finally remark that all these results do not 
depend sensitively on the details of the calculation, as is best seen by comparing the 
massive gluon calculations with those based on quasi-periodic boundary conditions. 

As mentioned at the outset, the specific choice of the Savvidy magnetic field 
in the gauge-fixed formulation we have pursued in this work has quite a different 
meaning to that in background field calculations. Because the choice of gauge has 
fixed both the choice of the diagonal color direction and the orientation of the z-axis 
in space, each possible choice of the colour and spatial orientation of the external 
field represents physically different cases. A natural extension of this work then is 
to perform similar calculations for such different orientations of the chromomagnetic 
field, both in internal and in Lorentz space. This program is nothing less than the 
mapping out of the physical directions of the effective potential. Similarly, it would be 
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interesting to study, in this gauge-fixed formalism, radically different external fields, 
such as self-dual fields [ 18| , which would include both magnetic and electric fields, or 
stochastic |19[ fields. Work in these direction is in progress. 
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Appendix: Landau Levels in the Axial Gauge 



The relevant part of the effective Lagrangian is 



(A.4i; 



with the definitions (A a denotes the color 3 component, i.e. the background potential) 



igf 



a/3 



d a + igA c 
d a ,d p 
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(A.42) 



Greek indices run from to 2. The Landau levels are most simply obtained from the 
classical Euler-Lagrange equations, 



(dpdP + d 3 d 3 + M 2 ) $ a - d a d^ + 2 [d a , dp] ^ = . 



(A.43) 



Applying d a to Eq. ( |A.43| ) and specializing to a constant magnetic background field, 
we can derive the constraint 

d a $ a = (A.44) 
which simplifies the Euler-Lagrange equations to 



(dpdP + d 3 d 3 + M 2 ) $ a + 2 [d a , dp] & = 



(A.45) 



Interestingly, Eq. (|A.44 ) which is usually imposed as background gauge condition 
can here be shown to be a consequence of the equations of motion. This is the reason 
why we recover eventually the standard results, although the background field in the 
gauge fixed framework has a somewhat different meaning. Let us describe the color 
neutral, constant magnetic field in 3-direction by the potential 



A 1 = -x 2 B 



A° = A 2 







and assume B > from now on. Then, the ansatz 

$ a ( X ) = e -i(^°-P^-^ 3 )^( x 2^ 



(A.46) 



(A.47) 
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converts Eq. (|A.45[) into 



-J 2 + E% + K. . 

-uj 2 + E 2 k + K. . 2igB 

-2igB -J 2 + E 2 k + h Ko , 




with 



E k = Vk 2 + M 2 (A.49) 
and the harmonic oscillator Hamiltonian characteristic for Landau levels, 



h h . . = -df + g 2 B 2 e, ^ = x2 + Jb~- (A ' 50) 

The constraint now assumes the form 

-itutp + igB£<p 1 + %? 2 = . (A.51) 

Introducing normalized eigenfunctions of /ih.o. via 

h hM .R n (Z)=gB(2n+l)R n (£)=n n R n (£) n = 0,1,2,... (A.52) 

and using the standard algebraic treatment of the harmonic oscillator to deal with 
the constraint, it is then easy to identify the following solutions of Eqs. (|A.48|) and 

<p{g)=M\ 1 Uo(£), uj 2 = E 2 k -gB 



i ^2g~B^+lR n (i)/uj 

Rn + i{Z) |, uj 2 = E 2 k + gB(2n + l) , n = 0,l,2,... 



V2g~B^iR n (0/uj 

<p(t)=M\ Rn-i(0 I, uj 2 = E 2 k +gB(2n+l) , n = l,2, 



(A.53) 



-iRn-l{t) 

The wavefunctions are given for the sake of completeness only. The eigenvalues uj are 
infinitely degenerate with respect to p; they enter into Eq. (|TTD of section 3, together 
with the appropriate boundary condition. 
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Figure Captions 



Fig. 1 Energy density of SU(2) Yang-Mills vacuum vs. magnetic field, for 3 different 
extension parameters L, in units where Ams = 1- Full curves: real part of 
renormalized energy density, Eq. (pT8|). Circles: quadratic approximation, Eq. 

Fig. 2 Regions of stability and instability in the (L, x) plane, where x is the angle for 
quasi-periodic boundary conditions. To the right of the circles, thermodynamic 
instability; above the solid line, magnetic instability. 

Fig. 3 Energy density and pressure, normalized to Stefan-Boltzmann values, vs. tem- 
perature in units of Ams- The curves correspond to the minimal model with 
quasi-periodic (massless) gluons, and should be regarded as upper bounds at 
high T. 

Fig. 4 Interaction measure e — 3P vs. temperature for quasi-periodic (x) and massive 
(M) gluons, compared with lattice calculation of Ref. ||. Units have been 
converted to GeV. 

Fig. 5 Regions of magnetic stability and instability in the (L, M) plane, where M is 
the gluon mass for massive (periodic) gluons. Solid curve: solution of Eq. fl3~2"|) . 
To the left of the circles, the system is unstable at all L. Units: A MS = 1. 

Fig. 6 Like Fig. 3, but for the minimal model with massive (periodic) gluons. 

Fig. 7 Electric screening mass of gluons as function of temperature, in units of Ams- 
Dashed curve: standard perturbative result with one loop running coupling 
constant, M: massive gluons, x : quasi-periodic gluons. 
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